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Abstract 

Aim  of  this  work  is  to  focus  attention  on  how  the  concentration  polarization,  in  all  electrochemistry  applications,  is  often  formulated  on  the 
basis  of  only  thermodynamic  statements  and  how  it  should  be  written  in  a  kinetic  more  consistent  expression.  Moreover,  the  application  of  this 
study  to  the  analysis  of  molten  carbonate  fuel  cell  (MCFC)  behaviour  has  been  discussed  with  reference  to  diffusion  limited  operating 
conditions. 
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1.  Introduction 

The  study  of  polarisation  losses  in  electrode  kinetics  is 
one  of  the  main  subjects  of  electrochemistry  and  of  major 
interest  today  considering  the  ever-increasing  worldwide 
interest  in  many  practical  applications  of  electrochemistry, 
such  as  fuel  cells. 

The  irreversibilities  occurring  in  an  electrode  reaction  are 
usually  summarised  in  three  polarisation  losses  to  be  deducted 
from  the  ideal  Nemst  potential:  activation  polarization,  ohmic 
polarization,  and  concentration  polarization  [1]. 

Activation  and  ohmic  polarisation  can  be  expressed  by 
theoretical,  semi-empirical  or  empirical  equations  which 
depend  on  materials,  temperature,  pressure  and  gas  compo¬ 
sition  [2]. 

Concentration  polarisation  is  very  important  when  the 
electrode  is  working  at  high  current  densities  and  it  is  related 
to  transport  phenomena  of  each  gaseous  reactant  in  the 
complex  structure  of  the  electrodes  (diffusion,  solution, 
adsorption,  etc.).  The  literature  usually  explains  this  voltage 
decay  by  means  of  thermodynamic  considerations  [1,3,4]. 

However,  it  is  worth  recalling  that  only  ideal  or  maximum 
efficiency  depends  on  electrochemical  thermodynamics, 
whereas  the  effective  efficiency  depends  on  electrode 
kinetics. 
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By  correctly  taking  thermodynamic  constraints  into 
account  in  a  kinetic  analysis  in  this  work,  we  wish  to  put 
in  evidence  the  connection  between  a  consistent  concentra¬ 
tion  polarisation  expression  and  the  usual  thermodynamic 
formulation. 


2.  Concentration  polarisation  on  the  basis  of 
thermodynamic  statements 

In  the  following  we  introduce,  with  a  formulation  of  wide 
validity,  the  thermodynamic  concentration  polarisation 
(AEtc)  usually  used  without  details  in  literature  [1,3,4]. 

We  will  refer  to  a  generic  heterogeneous  electrochemical 
reaction 

^  VjMf  +  vee~  =  0  (1) 

where  the  species  M,  migrates  to  or  from  the  electrodes. 

On  the  bases  of  the  difference  of  Nernst  voltage  when  the 
electrochemical  reaction  occurs  (that  is  at  electrode  surface 
concentration  cis)  and  when  no  current  is  supplied  (that  is  at 
bulk  concentration  cib),  the  thermodynamic  voltage  decay 
due  to  reactant  transport  phenomena  can  be  written  as 
follows: 

A Etc  =  Es  -  Eb  =  V  v,  In  ( — )  (2) 

veF  ^  \cibj 
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Nomenclature 

Ci 

reactant  concentration  of  the  component 
i  (mol  m~3) 

e- 

electron 

E 

voltage  (V) 

AE 

polarisation  (V) 

F 

Faraday’s  constant  (A  s  mol-1) 

J 

current  density  (A  m-2) 

k!  and  k" 

kinetic  constants  (A  m  mol-1) 

Kc 

mass  transport  coefficient  (m  s-1) 

Ke 

thermodynamic  equilibrium  constant 

Mt 

chemical  component 

R 

gas  constant  (J  mol-1  K-1) 

s 

real  positive  number  in  Eq.  (14) 

T 

temperature  (K) 

Zi 

charge  of  the  component  i 

Greek  letters 

ae  and  /)e 

reaction  rate  orders  related  to  e-  transfer 

coefficients 

a/  and  /?,■ 

reaction  rate  orders  related  to  the  component  i 

7 

thermodynamic  consistency  parameter  in 

Eq.  (14) 

n 

dimensionless  polarisation 

Ve 

stoichiometric  coefficient  of  e- 

Vi 

stoichiometric  coefficient  of  component  i 

Subscripts 

an 

anodic 

act 

related  to  reaction  activation 

b 

bulk 

cat 

cathodic 

C 

related  to  concentration  polarisation 

S 

electrode  surface 

0 

reference 

L 

limiting  value 

T 

thermodynamic 

ohm 

related  to  ohmic  resistance 

TOT 

total 

Moreover,  the  fraction  cis/c,- b  can  be  expressed  by  means  of  a 
variable  with  more  direct  physical  meaning  and  easier  experi¬ 
mental  evaluation:  the  local  limit  current  density  JL,  that  is  the 
maximum  current  reachable  when  the  limiting  reactant  con¬ 
centration  on  the  electrode  surface  approaches  zero. 

In  fact,  on  the  basis  of  Faraday  and  Fick’s  laws,  the 
following  equation  can  be  derived  for  current  density  [5]: 

J  =  FKaf  (Cft-Cfc)  (3) 

where  KCi  is  the  mass  transport  coefficient. 

The  limiting  current  density  J L  imposed  by  reactant  i 
corresponds  to  c,s  =  0  and  so  it  is 

Ju  =  FKCl-clh  (4) 

V; 


Inserting  Eq.  (4)  in  Eq.  (3),  it  is  possible  to  determine  c,s/clh 
as  a  function  of  JLi 


Cis  _  j  J 

C/h  J  Li 


(5) 


and  so  the  concentration  polarisation  can  be  written  as 


.  FT  v-  , 

AEtc  =  —  >  v;  In 

Ve 


(6) 


Moreover,  the  limiting  effects  of  only  one  reactant  (i  =  L) 
will  be  significant  in  comparison  to  the  other  ones,  so,  in 
dimensionless  form,  we  can  write 


rhc  ~ 


AEtcF  ^  vL 
RT  ~  ve 


In 


(7) 


with  JL  equal  to  the  minimum  among  the  JLi. 

However,  this  approach  is  not  a  completely  coherent 
combination  of  thermodynamic  and  kinetic  statements.  In 
fact,  A ETc  is  a  purely  thermodynamic  parameter  here,  just 
like  Eb,  whereas  we  should  be  interested  in  a  non-equili¬ 
brium  configuration  at  /  ^  0. 


3.  Concentration  polarisation  on  the  basis  of  kinetic 
statements 


The  Butler- Volmer  expression  for  current  density,  usually 
used  in  studying  electrochemical  reactions  [6-8],  is  here 
used  as  starting  point  in  order  to  deduce  the  concentration 
polarisation  directly  from  kinetic  statements  and  not  only 
thermodynamic  ones.  So,  we  assume  [9] 


J  Job 


(8) 


where  the  exchange  current  JQb  is 

Job  =  k'\{cl^FlRT  =  k"n4e-^/sr 

i  i 


(9) 


and  assuming  t]  =  (Es  —  Eb)F/RT. 

At  /  =  0  the  following  thermodynamic  consistency  con¬ 
straints  can  be  obtained 


Pi  ~  «/  =  &  +  °=e 

Vi  ve 


(10) 


k' 

k" 


k:  =  II 


vts  -vesEbF/RT 

Cib  e 


(ID 


where  s  is  a  real  positive  number.  Eqs.  (10)  and  (11)  are  a 
straightforward  generalisation  of  electrochemical  reactions 
of  similar  statements  of  common  use  in  chemical  kinetics, 
where  the  difference  in  the  sign  for  ae  in  Eq.  (10)  is  due  to  the 
opposite  direction  of  the  chemical  specie  e”  and  conven¬ 
tional  electrical  current. 

Moreover,  Ji)b  is  dependent  on  the  bulk  concentration  clb 
according  to  Eqs.  (9)  and  (11) 

Job  =  k*l[cl 


(12) 
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with 

fc*  —  idk~aelv  o 

(13) 

and 

v,ae 

Ji  =  a.'  + - 

(14) 

The  consistency  Eqs.  ( 1 0)— ( 1 4)  highlight  the  number  of 
degrees  of  freedom  of  the  kinetic  parameters;  when  the 
dependence  of  the  current  on  concentrations  is  known 
(namely  k*  and  yt),  together  with  stoichiometry,  the  knowl¬ 
edge  of  other  two  parameters  is  sufficient  for  the  determina¬ 
tion  of  the  number  s  and  of  all  the  orders  of  reaction  left. 

In  particular,  it  is  worth  noting  that  the  usual  position 
s  =  1  is  consistent  with  only  one  possible  stoichiometric 
choice. 

Now,  taking  into  account  Eqs.  (6),  (10)  and  (12)  the 
electrochemical  kinetic  can  be  expressed  as 

(>-£)> 

(15) 


i-n  i 


j 

Ju 


A  comparison  between  kinetic  dimensionless  polarisation 
0/)  from  Eq.  (15)  and  the  thermodynamic  concentration 
polarisation  (f]TC)  as  defined  by  Eq.  (7)  can  be  obtained  by 
considering  only  the  actual  contribution  of  limiting  reactant 
L  and  then  eliminating  the  current  density  J. 

In  particular,  after  some  manipulation  we  obtain 

JL(]  _  Q',e'hc/VL^  g-Vet/rcVi/Vz.  _  _  g- V„s(>/-t;rc)  j 

(16) 


with  JL  equal  to  the  minimum  among  the  JLi. 

Studying  this  equation,  we  can  observe  that  rj,  involving 
all  the  irreversibilities,  is  higher  than  i]TC,  which  is  related  to 
the  only  reactant  diffusion. 

In  particular,  the  linear  contribution  when  t]  —  0  is 

,i<o  <n) 


More  interesting  is  the  value  of  polarisation  when  diffusion 
phenomena  are  rate  determining,  that  is  when  the  electrode 
is  working  near  limiting  conditions  with ;/  — >  oo.  In  this  case 
the  value  of  i]  identifies  the  concentration  polarisation  tjc 
consistent  with  kinetic  statements  and  is 


ii  =  nr  =  —  In 


+  Itc 


0LevL)  aevt 

(18) 


so  a  difference  between  tjc  and  ;;rc  remains  when  yL  ^  0, 
that  is  uL/ae  /  -vL/ve. 

In  the  light  of  these  considerations,  the  concentration 
polarisation  has  to  be  written 


A  Ec  =* 


cllRT 
- In 


«e  F 


(19) 


and  is  a  function  not  only  of  JIJL,  but  also  of  the  coefficients 
ae  and  aL,  which  are  not  known  a  priori  like  the  stoichio¬ 
metric  coefficients  ve  and  iy. 

Fortunately,  on  the  condition  that  the  limiting  current  is 
well  determined,  the  difference  between  thermodynamic 
and  kinetic  formulations 

1c  ~  Itc  =  A Ec  -  A ETC  =  _  olA  _  l  =  _  (2Q) 

Itc  A  Etc  aev L  aev l 

is  not  very  important.  Also  when  yL  is  significantly  different 
from  zero,  the  consistent  kinetic  polarisation  and  thermo¬ 
dynamic  one  show  differences  in  the  VI J  performance  curve 
only  concerning  the  shape  of  the  “knee”  in  the  transition 
from  linear  behaviour  to  limit  current. 


4.  Application  to  molten  carbonate  fuel  cells 


The  definition  of  fuel  cell  performance  is  a  key  point  in  the 
development  of  these  electrochemical  devices  and,  even  if 
fuel  cells  are  usually  said  to  be  working  at  current  density 
much  lower  than  the  limiting  one,  the  operative  purpose  of 
maximising  the  fuel  utilisation  factor  approaches  the  diffu¬ 
sion  region,  so  that  the  cells  work  under  apparently  safe 
conditions,  but  with  possible  local  running  in  diffusion  limit 
conditions  [10,11]. 

This  has  been  shown  for  example  in  the  case  of  molten 
carbonate  fuel  cells  (MCFCs),  demonstrating  the  impor¬ 
tance  of  a  correct  evaluation  of  concentration  polarisation, 
that  we  will  discuss  below. 

In  MCFCs  the  electrochemical  reactions  are  the  following: 
H2  +  CO32-  ->  H20  +  C02  +  2e”  anode  (2 1 ) 

C02  +  j02  +  2e~  — >  CO32-  cathode  (22) 

Gaseous  reactants  are  then  anodic  H2  and  cathodic  C02  and 
02.  Each  of  these  has  to  pass  through  the  entire  complex 
electrode  structure  to  reach  the  reaction  site  and  this  con¬ 
tributes  to  concentration  polarisation. 

In  Fig.  1  a  simplified  scheme  is  presented  for  a  reference 
planar  MCFC.  In  this  picture  main  gas  flow  and  local 
transport  phenomena  directions  are  shown:  concentration 
polarisation  effects  will  have  to  be  calculated  for  each 
particular  position  along  the  cell  channels. 

In  particular,  if  the  cell  voltage  V  is  locally  written  as 

V  =  EScx  —  ESan  =  (£bca,  +  AEactcat  +  AEohm^,  +  AEccat) 

_  (^ban  F  AEactan  +  A£0hman  +  A Eq^  (23) 


we  obtain  A ECtot  =  A £Ccat  -  ATran  with  A ECca  <  0  and 
AECaa  >  0. 

The  concentration  polarisation  can  be  written  simply  as  an 
algebraic  sum  of  terms 


RT 

ATc'tot  - 


F 


aco2 


In  1  - 


J 


In  1  - 


JLCh 


J  LCCh 
_  OH; 


In  1  - 


■>lh2/  . 


(24) 
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OXIDANT  (main  direction) 

Fig.  1.  Simplified  scheme  of  a  reference  planar  MCFC. 


where  Ji h2  is  related  to  anodic  reaction,  while  both  Jlco2  and 
Jlo2  to  the  cathodic  one,  so  that  the  higher  among  these  last 
will  result  negligible  as  a  function  of  the  operating  condi¬ 
tions. 

Now,  evaluation  of  A Ecmr  required  the  knowledge  of  the 
reaction  orders  as  well  as  the  aean  and  aecat  coefficients, 
while  all  these  data  are  neglected  in  the  purely  thermody¬ 
namic  analysis. 

However,  experimental  data  collected  from  the  authors 
on  single  and  stacked  cells  operating  in  the  pressure  range 
1-3.5  atm  show  that  anodic  H2  and  cathodic  C02  concen¬ 
trations  strongly  affect  thermodynamic  voltage  in  accor¬ 
dance  with  Nernst  prediction,  but  do  not  affect  performance 
curve  slope  [2,10].  So,  if  we  assume  empirical  kinetics 
independent  of  the  concentrations  of  H2  and  C02,  we  can 


write  yHo  =  yCOl  =  0.  This  statement  means  that  the  order  of 
reactions  can  be  assumed  equal  to  the  opposite  of  stoichio¬ 
metric  coefficients. 

On  the  contrary,  the  earlier  mentioned  tests  [2,10]  showed 
empirical  kinetics  depending  on  the  concentration  of  02 
according  to  an  exponent  y0l  =  0.67. Then,  remembering 
that  vean  =  2,  vecat  =  -2,  vq2  =  -0.5  and  vHl  =  vCo2  =  -1, 
we  can  write 


A^Ctot 


RT 

F 


x  In 


(25) 


So,  the  evaluation  of  concentration  polarisation  should  also 
need  the  knowledge  of  the  aecat  coefficient.  This  is  usually 
considered  to  be  equal  to  1,  on  the  basis  of  the  current 
assumptions  s  =  1  and  ae  =  fic  (some  authors  maintaining 
s  =  1  distinguish  ae  from  fle  [6]).  Nevertheless,  an  expres¬ 
sion  of  concentration  polarisation  different  from  the  thermo¬ 
dynamic  one  is  always  obtained. 

However,  generally  02  diffusion  effects  are  negligible: 
MCFCs  always  work  with  air  fed  in  large  excess  to  the 
cathode,  that  is  under  conditions  in  which  the  density  current 
J  is  much  lower  than  the  limiting  current  due  to  oxygen 
transport  phenomena.  So,  the  second  term  of  Eq.  (25)  can  be 
neglected  and  the  concentration  polarisation  deduced  from 
the  thermodynamic  demonstration  is  acceptable,  but  only  on 
the  basis  of  certain  conditions. 

On  the  basis  of  these  considerations,  the  following  con¬ 
centration  polarisation  equation: 


AE'ctot 


RT 

2F 


In 


+  In 


(26) 


seems  to  be  useful  for  MCFC  study. 


Fig.  2.  Comparison  between  MCFC  experimental  and  simulation  data. 
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So,  in  order  to  verify  this  assumption,  the  Eq.  (26)  has 
been  integrated  into  a  detailed  MCFC  model  [2]  which 
allows  the  evaluation  of  global  MCFC  performance  thanks 
to  local  balances  of  mass,  energy  and  momentum. 

Fig.  2  shows  an  example  of  the  very  satisfactory  agreement 
obtained  between  simulation  results  and  experimental  data 
collected  on  a  laboratory  single  cell  of  Ansaldo  fuel  cells 
(anode:  C02  1.4  x  10“5  mol/s,  H2  5.7  x  10~5  mol/s,  H20 

2.3  x  ICC5  mol/s,  N2  16.6  x  10“ 5  mol/s;  cathode:  C02 

14.3  x  1CT5  mol/s,  N2  125  x  1(T5  mol/s,  02  15.1  x  1(T5 
mol/s).  Other  examples  at  different  operating  conditions 
are  available  in  [11]. 

5.  Conclusions 

Literature  usually  presents  a  concentration  polarisation 
expression  derived  only  from  thermodynamic  considerations. 

In  this  paper,  we  have  proposed  a  more  general  expression 
taking  account  of  the  thermodynamic  constraints  in  a  kinetic 
analysis.  The  difference  between  the  two  approaches  can  be 
particularly  evident  when  limiting  operating  conditions  are 
considered  and  reaction  orders  are  not  equal  to  stoichio¬ 
metric  coefficients. 

On  the  other  hand,  we  noted  that  thermodynamic  results 
can  be  acceptable  in  MCFC  application  thanks  to  a  parti¬ 
cular  combination  of  conditions. 
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